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The statistical properties of the intermittent signal generated by a recent model for self-organized- 
criticality (SOC) are examined. A successful comparison is made with previously published results 
of the equivalent quantities measured in the electrostatic turbulence at the edge of a fusion plasma. 
This result re-establishes SOC as a potential paradigm for transport in magnetic fusion devices, 
overriding shortcomings pointed out in earlier works [E. Spada, et al , Phys. Rev. Lett. 86, 3032 
(2001); V. Antoni, et al , Phys. Rev. Lett. 87, 045001 (2001)]. 
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Self-organized criticality (SOC) is believed to account 
for the behaviour of several extended non-equilibrium 
systems exhibiting bursty activity with long-range cor- 
relations in space and time 0. In the past years, SOC 
was candidated as a paradigm for the understanding of 
anomalous transport of energy and matter in magneti- 
cally confined fusion plasmas j2J. Within the huge com- 
plexity of transport in fusion devices, in fact, a num- 
ber of features were identified, which could easily be 
cast into the framework of SOC systems: from the ex- 
istence of critical average gradients and profile resilience, 
to the power-law power spectra of plasma parameters 
fluctuations (density, temperature, magnetic field, ...). 
Many properties of SOC numerical models are satisfac- 
torily compared against experimental data [E 01 IE IE : 
e.g. transport barriers have been recently reproduced us- 
ing SOC models 0; non- locality, which is intrinsic in 
SOC, appears to be a possible ingredient for core trans- 
port ; numerical simulations of edge turbulence gave a 
phenomenology SOC-like 0. 

However, recently further analysis showed some irre- 
ducible discrepancies between numerical and experimen- 
tal time series |2l llOj : these claims have by now con- 
firmed in almost all devices. Essentially, numerical time 
series (normally coming from the Bak-Tang-Wiesenfeld 
(BTW) model |l| , the sandpile d la Kadanoff [ll| , or the 
Hwa-Kardar "running" sandpile were not found to 
display realistic temporal correlations. It was pointed 
out that experimental data display distribution of wait- 
ing times with power-law tails, and that intermittency 
appears as an ubiquitous property of plasma time series: 
a behaviour not registered in the analyzed SOC mod- 
els. The first paper addressing the former point in the 
context of laboratory plasma was Ref. 0, carrying on an 
analysis of the probability distribution function (PDF) of 
waiting times between bursts in density fluctuations mea- 
sured on the RFX Reversed Field Pinch Experiment . 
Experimentally, it was found a power-law curve with an 
exponent about —2, while the considered SOC models 
had exponential-like PDFs. Later, it has become clear 



that the statistics of waiting times is not truly a strin- 
gent test about the existence of SOC, because several 
SOC models have also non-exponential PDFs of waiting 
times (e.g., see 0,0] and references therein). 

Intermittency in turbulence implies the lack of self- 
similarity between time scales It has been charac- 
terized in plasma experiments by looking at the PDF of 
signal differences (or with the most sophisticated tool of 
the continuous wavelet transform 0]) at different time 
scales 0, . It turns out that the shape of the PDFs 
does not collapse to a single curve, irrespective of the 
time scale. It was found the presence of non-Gaussian 
PDFs in the RFX signal, approximately stretched expo- 
nentials, at the smaller r's (higher frequency), gradually 
recovering the Gaussian shape as r increased. Again, this 
is not the case in the analyzed SOC models. 

SOC alone is a paradigm, which needs to be imple- 
mented into specific models in order to provide verifiable 
predictions. Part of the extracted predictions, hence, 
will be intrinsic to SOC, while others will depend on 
the model. Normally, as long as a specific model fails 
to account for some empirical evidence, one may devise a 
variant of it that is able to cure the specific shortcoming, 
but the improved accuracy of the final model is obtained 
at the expenses of some loss of generality. In our case, the 
simplicity of the SOC paradigm would be spoiled and ob- 
scured by model-specific features. A good option would 
be one where the added features reflect a property that 
the true physical system is anyway likely to possess. 

Previously, turbulent features were already found in 
time series of "waves" in the BTW model Oil . The same 
features were later found in experimental time series of 
solar flares (2jj. Waves are decompositions of avalaches, 
and by definition, within an avalanche, they start from 
the same site. Thus, a message from this fact is that 
spatial features of the drive are important, affecting the 
temporal correlations in the activity. 

Sanchez, Newman and Carreras advanced several sug- 
gestions, to make compatible the findings of Refs. [£| |lfj 
within a SOC-like scenario. The first proposal |2l| was 
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FIG. 1: A sample of the time trace of the output of the model, 
showing its bursty behaviour. 



FIG. 2: Power spectrum of the sinthetic signal for both pe- 
riodic and open-closed boundary conditions. Overplotted for 
comparison is a 1// curve. 



that of a mixed avalanche-diffusive transport, where the 
self-similarity of timescales would be naturally broken by 
the second contribution. A later suggestion 22] hinted 
to the possibility of inducing modifications to the wait- 
ing times distribution through correlations built into the 
drive, hence remaining within a pure SOC framework. 
However, the authors pointed out that introducing this 
sort of explicit correlations into a SOC model could be 
not trivial, since the generating mechanism is likely to 
be specific to any physical system, and therefore scarcely 
useful. We show that indeed it turns out to be a pes- 
simistic view: a very simple and general correlation 
mechanism can still generate a very complex time series 
and output. 

We are going to present in this work a SOC model that 
is able to reproduce the findings of papers 0, El • This is 
obtained by postulating the existence of a correlation be- 
tween successive locations of the fueling. This requisite 
appears rather plausible in a model that attempts to sim- 
ulate the physics in a plasma device: particle and heat 
sources in these environments are more or less localized, 
in any case not uniformly randomly distributed. 

The SOC system is visualized as a standard 1-dim lat- 
tice automaton of length L, very similar to the sand- 
pile studied by Kadanoff et al The choice of the 
model and of the boundary conditions will allow a cor- 
respondence between its profile h and ordinary profiles 
in laboratory plasmas. Each site i (1 < i < L) of 
the lattice holds hi "grains". Different choices for the 
boundary conditions can be taken: in |l4| periodic ones 
(/iL+i = hi, ha = /il) were used. In this work, we will 
use one open boundary (/il+i = 0), while the other one 
is closed (ho — hi). Stable pairs of adjacent sites 
fulfill the local stability condition if \hj — hi\ < H, where 
the constant H is a threshold. A local instability is re- 
solved by a toppling, which consists in moving a grains 
from the upper to the lower site. This lower site in turn 
can become unstable, and chain reactions of transport are 



possible (avalanches). All instabilities in the systems are 
updated in parallel, and this procedure is iterated until 
the avalanche ends because all sites returned stable. The 
whole process takes place in one "time step". Instabili- 
ties arise because the system is driven out of equilibrium: 
a new grain is added at each time step at a position i: 
hi — ► hi + 1. Correlations are introduced at the stage 
of choosing a new site i' for deposition at time step t, 
provided that at the previous step (t — 1) the grain had 
been deposited on site i. In earlier models \K HH fl^| . 
the choice is completely random and uncorrelated from 
the value i. Here, on the contrary, we make the opposite 
hypothesis of a strongly correlated diffusing dynamics: 
with equal probability, i' = i + 1 or i' — i — 1. Reflecting 
boundary conditions for the drive are used to avoid i' to 
fall outside the lattice. 

We do not expect qualitatively different results by 
varying slope-related quantities like H and a. Hence, 
we will keep them fixed from now on: H = 4, a = 2. We 
shall see, instead, that some results are not L-invariant. 
For reference, we will use L = 512. In what follows we 
will monitor for diagnostic purposes the total activity s, 
i.e. the total number of topplings within each avalanche. 
This activity is naturally identifiable with the quantity of 
matter or heat delivered to the device's walls, hence close 
to experimental investigation. All the simulations were 
carried out by starting with an empty system and let- 
ting it to load until a stationary critical phase is reached. 
Then, recording phases of several 10 6 time steps were 
performed. One sample of the typical time series is given 
in Fig.HJ 

The power spectrum S(f) is a statistical measure that 
can be easily computed from data sets, and gives funda- 
mental information about time correlations built into the 
signal. Hence, an important issue is that about its scaling 
properties. The signature of SOC models is power-law 
spectrum. Most experiments (including also "numerical 
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experiments" ) claim to find a power-law, or possibly more 
than one over different frequency ranges, with a variance 
of the exponent. Its absolute value, in several works, is 
higher than unity over most of the sampled frequency 
range, usually ranging between 2 and 3 |23j. There are 
however some exceptions 0, 0, 0, EEj where a 1// scal- 
ing appears to be recovered. It is possible that in the 
first set of papers the 1// region exists too, but its width 
is too small to be detectable, and what was measured 
was instead just the high-frequenc y pa rt of the spectrum 
(see, in this regard, Fig. 3 of Ref. |25j|). 

In Fig. |21 we show the power spectrum from the model, 
for two choices of boundary conditions (we study both 
because the choice of boundary conditions does affect the 
spectrum). In the periodic boundaries case, a 1// slope 
is clearly discernible. In the other case, the trend still 
appears but lesser clear. This is quite remarkable a result: 
it is well known that the power spectrum of a random- 
walk-like signal is a l// 2 curve. Hence, the random- 
walk dynamics of the driver couples to the rules for the 
stability of the sandpile to produce an highly nontrivial 
output. 

The first quantity of interest is the PDF of waiting 
times between bursts, P(t w ), shown in Fig.0 We define 
as bursts the points s(t) for which s(t) > 3 x (s), with 
(s) the average value over the sample. Most of the P(t w ) 
curve is accurately fitted by a power-law with exponent s=s 
—2; interestingly, it is quite close to experimental values 
from RFX. The appearance of a non-Poissonian PDF is 
related to the existence of time correlations inherited by 
the driving [T^ |. 

So far, we have shown that the present model fulfills 
the questions raised in work 9] as long as the scaling 
of waiting times is concerned. The issue of the depar- 
ture from self-similarity is positively addressed in Fig. 0] 
There, we plot the PDF of the wavelet transform of the 
signal at three time scales. The PDF has fat tails at the 
smaller time scales, gradually approaching a Gaussian 



FIG. 4: Distributions of coefficient of wavelets transform at 
three time scales. The heights of the PDFs have been scaled 
to their maximum value, and the widths to the standard devi- 
ation. The solid curve is a Gaussian with the same amplitude 
and variance, plotted for reference. Upper plot, the lattice 
length is L — 512; lower plot, for comparison, L = 128. 

shape at the longest time scales. The explanation we give 
for the breaking of the self-similarity is this: the seeding 
has an effective " diffusivity" D = (length step) 2 /(time 
step) = 1. After a time of order T = L 2 /D = L 2 the 
seeding has sampled the whole lattice, irrespective of the 
starting position. Thus, T plays the role of a correlation 
time: the system keeps some memory of its past for as 
long as this time. Semi-quantitatively, some confirma- 
tions of this statement may be found: I) in figure 0] we 
show that the convergence to Gaussian shape is faster 
with decreasing L. II) By sampling the system at a fre- 
quency / < 1/T one is actually sampling from indepen- 
dent realizations of the same system (Gibbsian ensem- 
ble). Hence, the scalar signal picked out must be equiv- 
alent to a random number from within a uniform finite 
distribution, whose power spectrum is white noise, / . 
Indeed, the power spectrum in fig. |21 begins to flatten 
towards an /° spectrum for / lesser than about 10~ 5 , 
consistent with the 1/L 2 w 4 x 10 -6 estimate. 

Although successful when benchmarked against all the 
sought statistical tests, the model here proposed is still 
quite minimal, and add-ons may be envisaged that can 
enrich the built-in physics without appreciably spoiling 
its simplicity. An option we are pursuing is that of im- 
plementing into the model the running sandpile version 
of Hwa and Kardar [12j . Indeed, a time scale common to 
both the drive and the system dynamics appears a suit- 
able feature for models of turbulent phenomena 

Nowadays, the intermittent character of turbulence at 
the edge of fusion devices is attributed to the existence of 
long-lived coherent structures moving on top of a back- 
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ground plasma (dubbed blobs), which are responsible 
for most of the transport to the wall H||. The precise 
mechanism of blob generation is self-regulating and in- 
volves the triggering of an interchange instability driven 
by pressure gradients with formation of a radially elon- 
gated structure, that ultimately is separated from the 
core plasma by the differential rotation and expelled to- 
wards the edge, with an accompanying relaxation of the 
mean profile 0, H3|. The affinities with SOC dynam- 
ics are apparent, and have been noticed by some authors 
(see, e.g. [11,113); however, to our knowledge, the under- 
standing of blobs' physics has not yet advanced enough 
to allow for quantitative comparisons with SOC models. 

Undoubdtedly, the overall edge transport comes from 
an intricate pattern where the physics of creation, de- 
struction, motion, and interaction of the coherent struc- 
tures between them and with the background, needs to 
be carefully taken into account in order to provide a cor- 
rect simulation of experimental evidence. It is there- 
fore still an open question whether the SOC paradigm 
would be comprehensive enough to account for all this 
physics. The evidence presented in this paper suggests, 
however, that further analysis is necessary before dis- 
carding SOC as a potential transport paradigm. Indeed, 
in the past there has been sometimes a tendency to see 
a dichotomy between SOC and magnetohydrodynamics 
(MHD) turbulence that, perhaps, is there not: attempts 
of developing cellular automata consistent with MHD 
and Maxwell's equations are actually found in literature 
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